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Abstract 

Low-energy scattering of D* and Di meson are studied using quenched lattice QCD 
with improved lattice actions on anisotropic lattices. The calculation is performed 
within Liischer's finite-size formalism which establishes the relation between the 
scattering phase in the infinite volume and the exact energy level in the finite 
volume. The threshold scattering parameters, namely the scattering length oq and 
the effective range rg, for the s-wave scattering in = 0~ channel are extracted. 
After the chiral and continuum extrapolations, we obtain: ag = 2.52(47)fm and rg = 
0.7(l)fm where the errors are purely statistical. Based on these results, we discuss 
the possibility of a shallow bound state for the two charmed mesons within the 
non-relativistic potential scattering model. It is argued that, albeit the interaction 
between the two charmed mesons being attractive, it is unlikely that they can 
form a shallow bound state in this channel. This calculation provides some useful 
information on the nature of the newly discovered resonance- like structure Z+(4430) 
by the Belle Collaboration. 
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1 Introduction 

Recently, a charged resonance-like structure ^+(4430) has been observed at 
Belle in the vr^' invariant mass spectrum of 5 — Kir^ip' decays [1]. This 
discovery has triggered many theoretical investigations on the nature of this 
structure [2|5|HI5|5f7f5f51[TU] . Since the invariant mass of the resonance is very 
close to the D*Di threshold, one possible interpretation is a molecular bound 
state formed by the D* and Di mesons [3|H] . To further investigate this pos- 
sibility, the interaction between D* and Di mesons becomes crucial. As is 
known, the interaction of two hadrons can be studied via the scattering pro- 
cess of the hadrons. Since the energy being considered here is very close to 
the threshold of the D*Di system, only threshold scattering parameters, i.e. 
scattering length and effective range ro, are relevant for this study. In phe- 
nomenological studies, the interaction between the mesons can be computed 
by assuming meson exchanges models. However, since the interaction between 
the mesons at low-energies is non-perturbative in nature, it is tempting to 
study the problem using a genuine non-perturbative method like lattice QCD. 

In this paper, we study the scattering threshold parameters of D* — Di sys- 
tem using quenched lattice QCD within the so-called Liischer's formalism, 
a finite-size technique developed to study scattering processes in a finite vol- 
ume pT|12fl3lll4fl5] . Within this approach, it is also feasible to investigate the 
possible bound state of the two mesons [T^IITB] . We have used improved gauge 
and fermion lattice actions on anisotropic lattices. The usage of anisotropic 
lattices with asymmetric volumes has enhanced our resolution in energy and 
the momentum. The computation is carried out in all possible angular mo- 
mentum channels, although only the = 0~ channel yields definite results. 
We find that, in this particular channel, the interaction between a D* and a Di 
meson is attractive in nature. The scattering length after continuum and chi- 
ral extrapolation is oq = 2.52(47)fm while the effective range is ro = 0.7(l)fm. 
Possibility of a bound state can also be addressed within Liischer's formal- 
ism. Our simulation results indicate that the two-particle system of the two 
charmed mesons resembles more like an ordinary scattering state rather than 
a shallow bound state. 



^ Work supported in part by NSFC under grant No. 10835002, No. 10675005 and 
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This paper is organized as follows. In Section [2], we briefly introduce Liischer's 
formalism and its extensions to the asymmetric volumes. In Section [3l we dis- 
cuss possible one-particle and two-particle interpolating operators and their 
correlation matrices are defined. In section HI simulation details are given and 
the results for the single- and two-meson systems are analyzed. After verify- 
ing the single-particle states, we extract the exact energy of the two-particle 
system. By applying Liischer's formula the scattering phases are extracted for 
various lattice momenta. When fitted to the known low-energy behavior, the 
threshold parameters of the scattering system, i.e. the scattering length 
and the effective range vq are obtained in the s-channel. We also discuss vari- 
ous interpolation and extrapolations which bring our results to the chiral and 
continuum limit. Based on our simulation results, the possibility of a bound 
state in this channel is discussed. In Section |5l we will conclude with some 
general remarks. 



2 Strategies for the computation 

2.1 Liischer's finite volume technique and its generalization 

Within Liischer's formalism, the exact energy eigenvalue of a two-particle sys- 
tem in a finite box of size L is related to the elastic scattering phase of the two 
particles in the infinite volume. Consider two interacting particles with mass 
mi and m2 enclosed in a cubic box of size L, with periodic boundary con- 
ditions applied in all three directions. The spatial momentum k is quantized 
according to: 



with n being a three-dimensional integer. Now consider the two-particle system 
in this finite box and let us take the center-of-mass frame of the system so that 
the two-particles having three-momentum k and — k respectively. The exact 
energy of the two-particle system in this finite volume is denoted as: Ei.2(k). 
We now define a variable k^ via: 



Note that due to interaction between the two particles, the value of k^ differs 
from its free counter-part k^ with k being quantizes according to Eq. ([1]). It 
is also convenient to further define a variable as: 






(2) 



= PLV(27r)2 . 



(3) 
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which differs from due to interaction. What Liischer's formula tells us is 
a direct relation of and the elastic scattering phase shift tan 5(g) in the 
infinite volume and it reads: [H] 



where Zoo{^', <f) is the zeta-function which can be evaluated numerically once 
its argument g^ is given. Therefore, if we could obtain the exact two-particle 
energy i?i.2(k) from numerical simulations, we could infer the elastic scattering 
phase shift by applying Liischer's formula given above. 

Here we would like to point out that, the above relation is in fact only valid 
under certain assumptions. For example, the size of the box cannot be too 
small. In particular, it has to be large enough to accommodate free single- 
particle states. Therefore, in a practical simulation, one should check whether 
this is indeed realized in the simulation. Polarization effects are also neglected 
which are suppressed exponentially by (9(e~'"^) where m being the single- 
particle mass gap. Also neglected are mixtures from higher angular momenta. 

In the case of attractive interaction, the lowest two-particle energy level might 
be lower than the threshold which then renders the quantity g^ < 0. The phase 
shift in the continuum, 5(g), is only defined for positive g^, i.e. energies above 
the threshold. When g^ < 0, it is related to the phase cr(g) via: 



where (— ig) > and the phase a{q) for pure imaginary g is obtained from 
5(g) by analytic continuation: tancT(g) = —i tan 5(g) piUITB] . The phase cr(g) 
for pure imaginary g is of physical significance since if there exists a true 
bound state at that particular energy, we have tan(T(g) = —1 in the infi- 
nite volume and continuum limit. In the finite volume, the relation above is 
modified as: [16] 



where the finite-volume corrections are assumed to be small. Therefore, for 
g^ < 0, we could compute tancr(g) from Monte Carlo simulations and check 
the possibility of a bound state at that energy. 

The above formulae apply to the case of a box with cubic symmetry. In 
real calculations, in order to have more accessible low-momentum modes, it 




(4) 




(5) 




(6) 
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is advantageous to use asymmetric volumes in the study of hadron scatter- 
ing [T7|18fl9j . If the rectangular box is of size L x (772-/^) x {rj^L), then Eq. (jlj) 
is modified to: 

tan6{q) = ^—— — -, (7) 



where the modified zeta-function Zoo(l, 1?^; ?72, ''^s) is the analogue of 

and its explicit definition can be found in Refs. [T7|IT8] . Similarly, for negative 

q^, the formula is modified to: 

tana(g) = -. (8) 



3 One- and two-particle operators and correlators 



Single-particle and two-particle energies are measured in Monte Carlo simula- 
tions using appropriate correlation functions. These correlation functions are 
constructed from corresponding interpolating operators with definite symme- 
tries. Since we are interested in the interaction between a D* and a Di meson, 
we need one-particle operators which would create a single D* and a single Di 
meson and two-particle operators which create both D* and Di from the QCD 
vacuum. Below we will first list these one-particle and two-particle operators 
and then proceed to discuss their correlation functions. 



3.1 One- and two-particle operators with definite symmetries 



Let us first construct the single meson operators for D*{2010)^ and Z)i(2420)° 
whose quantum numbers are 1~ and l"*", respectively. Just to simplify the 
notation, we will use Qi and Pi for these meson operators respectively, where 
i = 1, 2, 3 being the index to specify different spatial components. We use local 
interpolating fields as follows: 



Q,(x) = [(iVc](x), P,{x) = [cYl'u]ix) (9) 

where Qi{x) stands for L>*(2010)+ while P,{x) stands for Z)i(2420)°. A single- 
particle state with definite three-momentum k is represented by the Fourier 
transform of the above operators: 



Qi{t, k) = ^ Q,{t, x)e-^^ ■^ P,(t, k)=Y, P,{t, x)e-^'^-^. (10) 
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Obviously, the operators Qi{t, k) and Pi{t, k) fall into the vector representation 
of the rotational group S0{3) (i.e. their angular momentum quantum number 
is 1) in the continuum. 

On the lattice, the rotational symmetry group SO (3) is broken down to the 
corresponding point group. Usually, one utilizes an symmetric cubic box. In 
this case, the corresponding point group is the cubic group 0(Z). However, 
in order to access more non-degenerate low-momentum modes, it would be 
advantageous to use asymmetric box (although the lattice spacings in spatial 
directions are still symmetric). This is particularly useful for scattering pro- 
cesses, as advocated in Ref. [19]. Following this strategy, we have adopted a 
rectangular box of size L x {r]2L) x (rj^L) with t]2 = 1 and 773 7^ 1. In this case, 
the rotational group in the continuum is broken down to the basic point group 
D4. In what follows, we will construct operators that transform according to 
different irreducible representations (irreps) of the D4 group. 

The basic point group D4 has four one-dimensional irreducible representations: 
Ai, A2, Bi, B2 and one two-dimensional irreducible representation: E. With 
these notations, it is easy to verify that three components of an ordinary vector 
in the continuum, like Qj's and Pj's given above, now falls into two irreps: A2 
and E. In particular, we have the following decomposition rules: 

= Ai, 1 = E © A2, 2 = Ai®Bi®B2®E. (11) 

For the two-particle system formed by a D* and a Di meson, the quantum 
number of the two-particle system can be: = 0~,1~,2~. Now, we 
consider the vector space {Qi, Q2, Q3} ® {-Pi, P2, P3}, which is 9-dimensional. 
Using standard group-theoretical methods, it is easy to find out that this 9- 
dimensional vector space is made up of two copies of Ai, one copy of A2, Bi 
and B2 each and two copies of E. The basis operators of each irrep mentioned 
above are listed as follows: 

0(^0(1) (t) = Y^[Q^^t + l,-Ro k)Pi(t, R o k) 

ReG 

+ Q2it+l,~Rok)P2{t,Rok) 

+Qs{t + l,-Rok)P,{t,Rok)], 

0(^^)(^){t) = J2[Q,{t + l,-Ro k)Pi(t, R o k) 

ReG 

+Q2it + l,-Rok)P2it,Rok) 
-2Q3{t + l,-Rok)P3{t,Rok)], 

0(^^)(t) = J2 [Qiit + 1, o k)P2(t, R o k) 

ReG 

-Q2it + l,-Rok)Piit,Rok)], 
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Table 1 

The two-particle operators defined in Eq. (jl2p and their corresponding angular 
momentum quantum number J in the continuum. 





Two-particle operators 




C)(Ai)(i)(t) 


1" 




2 


o(^i)(2)(t),o(^i)(t),o(^2)(t),op(')(t),of)(')(t) 



[Ql(t + l,-i?ok)Pi(t,i?ok) 

R€G 

-Q2(t + l,-i?ok)P2(t,i?ok)], 

^ [gi(t + l,-i?0k)P2(t,i?0k) 

+(52(t+ l,-i?ok)Pi(t,Pok)], 

5] [gi(t + i,-Pok)P3(t,Pok) 

rgg 

-g3(t+l,-i?ok)Pi(t,Pok)], 

^ [g2(t + i,-i?ok)P3(t,Pok) 

-g3(t + i,-Pok)P2(t,Pok)], 
[gi(t + i,-Pok)P3(t,Pok) 

i?GG 

+g3(t + i,-Pok)Pi(t,Pok)], 
5] [g2(t + i,-Pok)P3(t,Pok) 

i?GG 

+g3(t + i,-Pok)P2(t,Pok)], (12) 

where k is a chosen three-momentum mode and G is the group and R & G 
is an element of the group. 0(^i)W(t) with z = 1,2 in this case designates 
different copies of the Ai representations occurring in the decomposition. Note 
that in the above definitions we have not included orbital angular momentum 
of the two-particles. Therefore we are only studying the s-wave scattering 
of the two mesons. This is sufficient for this particular case since near the 
threshold, the scattering is always dominated by s-wave contributions. Using 
the correspondence in Eq. (fTTI) . it is easy to figure out the continuum quantum 
numbers for these operators which are tabulated in table [H 

3.2 Correlation functions 

We then proceed to discuss one-particle and two-particle correlation functions, 
respectively. As already mentioned, in the lattice study of hadron-hadron scat- 
tering, one first have to make sure that asymptotically free one-particle states 



of = 
of = 
of )(^)(t) = 
of )(^)(t) = 
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are realized in the volume being considered. We therefore construct the one- 
particle correlation function C^{t, k) and C^{t, k) for the D*~^ and meson 
as: 



cQ(t,k) = (g,(t,k)g,(o,k)t) 

= -Ee-^'^-(7.75).,(7.75).. (xSr(-k)) (Xijir; 

X 

C^(t,k) = (P,(t,k)P,(0,k)t) 

= - E e-'^-hM7d.s (Xgr ^ (-k)) (Xiir y , (13) 

X 

where k is the three-momentum of a single meson. The quantities like X'^^^J^J^^^ 
stands for the quark propagator for a particular flavor /. For example: 



L J /3,a,x,t;7,6,y,0 

xiST^ =y:\M^'n'' . (14) 

y 

where we have assumed that the fermion matrix J\A^^^ satisfying: A^^-'^)^ = 
75A^(-^)75 for any flavor /. In the large temporal separation limit, the energy 
i?(k) of a single meson with definite three- momentum k can be extracted from 
the effective mass plateau of the corresponding correlation functions as usual. 

Next, we will discuss the more complicated two-particle correlation functions. 
Generally speaking, we need to evaluate a correlation matrix of the form: 

{o^PKt)opm, (15) 



where T and T' labels the irreducible representation of the group (i.e. T = Ai, 
A2, Bi, B2 and E for group ). However, as we show below, we do not need 
to calculate the whole matrix in Eq. (fT5l) . Since these point group represen- 
tations are all real, the hermitian conjugate of an operator transforms in the 
same manner as the original operator. Furthermore, since the QCD vacuum 
is invariant under any group transformations, it is therefore seen that, only 
the invariant sector (i.e. Ai sector), decomposed from the product of two ir- 
reducible representations: F (g) F', can make a non-vanishing contribution to 
the correlation matrix defined above. For the group D4, all irreducible repre- 
sentations are one-dimensional except E which is two-dimensional. Therefore, 
the direct products of two irreducible representations are particularly simple. 
For example, we easily verify that, the direct product of any two different 
one-dimensional irreducible representation cannot contain the Ai representa- 
tion while the direct product of any one-dimensional irreducible representation 
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with itself is exactly the Ai representation. It is also seen that, the direct prod- 
uct of any one-dimensional irreducible representation with E also contain no 
Ai components. We therefore only have to consider the combination E ® E 
which reads: 

E®E^Ai®A2®Bi®B2 (16) 

So two E operators can yield an invariant representation Ai. But the other 
three ingredients cannot contribute. From this discussion we conclude that, we 
only have to consider the case F = F' and in particular, iiV = E, we only have 
to consider something like: Z^Q=i(O^^^''"(t)O^^^(0)). However, one should keep 
in mind that, different momentum modes k do mix. This is what causes the 
scattering. The argument given above implies that, we only have to compute 
one correlation matrix in each channel. The size of this matrix is n x n where 
n is the number of momentum modes being considered. 

For the operator 0'''^^'>^^\t) , the correlation function is as follows: 

^il^^^^H^) = (oL^^^^^^^Woi^^^^^Ho)), (17) 

where m and n are indices for different momentum modes. We notice that this 
correlation function has a disadvantage in practical calculations. The sum- 
mation over the group element R in the definition of the operator 0^-^'^'^^^\t) 
cannot be absorbed into the source-setting when solving the propagators. This 
drawback can be cured by using a slightly modified operator: 

0(Ai)(l)(^)=^ E Q^(t+l,-i?'ok)P,(t,i?"ok) (18) 
i=l R',R"eG 

The difference of this operator as compared with the original operator is that, 
this operator contains also non-zero total three-momentum components. To 
be specific, those terms with R' ^ i?", will create states with non-zero total 
three-momentum. However, if we form the correlation function: 

^il^^^^n^) = (OL^^^^^^n^)Oi'^^^^^H0)), (19) 

then since the sink operator has total three-momentum zero, and the vac- 
uum also has total three- momentum zero, only the zero momentum terms 
in O^^i^*^^'(0) will contribute to the correlation function. That is to say, this 
will yield same correlation function as the original operator. However, using 
the operator O*^'^^'"-^-'(0) at the source has a big advantage. It will allow us to 
complete the summation over Rl and i?" in one step. As a result, instead of 
solving for the quark propagators for each i?, we only have to solve the quark 
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propagator once, with R being summed over and absorbed into the source 
definition. 



Implementing the trick mentioned above, the final result of correlation function 
according for the operator 0^^'^^^-^\t) is as follows: 



ReG i,j=i X 

■(E 4Sl'^(^''°q))-(^Sr'°^; 



(20) 



with m and n being momentum mode indices with corresponding three-momenta 
p and q, respectively; R E G being a group element of D^, X being the quark 
propagators with appropriate sources as defined in Eq. ([E 



Another important feature that has became clear from the above expression 
is that, the light quark propagators are needed for the zero momentum mode 
only. Different momentum modes enters the heavy quark propagators. Since in 
the quark propagator inversions, light quarks cost most of the computer time, 
this separation means that we only have to solve the most time-consuming 
part of the propagator, which is the light quark propagator, for vanishing 
three-momentum. Heavy quark propagators are needed for each momentum 
mode, however, it is not costly since the quark mass is heavy. 



Two tricks mentioned above, one being the reduction to non-degenerate mo- 
mentum modes, i.e. different three-momentum modes that are related by D4 
group transformations requires only one quark propagator inversion; the other 
being solving light quark propagators for zero momentum only, have offered 
us enormous amount of acceleration in the calculation. For example, taking 
highest three-momentum up to (1,1,0), the number of different three-momenta 
is 21. But if we only count the non-degenerate momentum modes, it is only 
6, gaining more than a factor of 3. Now that we only have to compute the 
zero momentum mode for the light quark, this gives again a factor of almost 
6 (neglecting computer time for heavy quark inversions). Altogether, we ex- 
pect a factor of about 15-20 gaining in the speed of the simulation. Note also 
that, these tricks are generally applicable for any type of calculations involving 
mesons with one heavy and one light quark. 
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4 Simulation details 



4-1 Lattice actions and simulation parameters 

The gauge action use in this study is the tadpole improved gauge action on 
anisotropic lattices: """""" 



S- 



i>j L 



5 TrP,;, 



1 TrRij 



1 TrRji 
~ 36 ~ 36 

4eTrPoi 1 ^TrRio' 



9 u2 



36 ut 



(21) 



where Pjj is the usual spatial plaquette variables and Rij is the 2x1 spatial 
Wilson loop on the lattice. The parameter Ug, which we take to be the 4-th 
root of the average spatial plaquette value, incorporates the so-called tadpole 
improvement and ^ designates the aspect ratio of the anisotropic lattice. The 
parameter /3 is related to the bare gauge coupling which controls the spatial 
lattice spacing in physical units. 

The fermion action used in this study is the tadpole improved clover Wilson 



action on anisotropic lattice whose fermion matrix is [23.24 J: A4xy 
with: 



xy 



1 3 

Axy = Sxy[- \- PtT (^oiFoi + PsicTuFu + (T23-F23 + 0-31-^31)] 

^^max i—i 

- E ^^[(1 - l,)U,ix)Sx+,,y + (1 + l.Wlix - i,)5,.,J (22) 
where the coefficients are given by: 



Vi 
pt 



2^' - 2' 



1 

2k 



2k r. 



Ps 



V 

2^' 



(23) 



Quenched gauge field configurations are generated using the conventional 
Cabbibo-Mariani pseudo-heat bath algorithm with over-relaxation. Quark prop- 
agators are obtained using the so-called Multi-mass Minimal Residual (M^R) 
algorithm, which can yield the propagators with different quark masses at one 
inversion [25]. Dirichlet boundary conditions are used in the temporal direc- 
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Table 2 

Simulation parameters in this study. All lattices have the same aspect ratio: ^ = 5. 





p = 2.5 


p = 2.8 


p = 3.2 


-^conf 


700 


500 


200 




0.4236 


0.4630 


0.50679 




0.732 


0.79 


0.89 




0.9305 


0.96 


1.0 


asifm) 


0.2037 


0.1432 


0.0946 


lattice 


8 X 8 X 12 X 40 


12 X 12 X 20 X 64 


16 X 16 X 24 X 80 




0.0577 


0.0598 


0.0595 


i^ud 
"'max 


0.0613 


0.0611 


0.0606 



tion for the fermion fields. Error estimates are made using the conventional 
jack-knife method for all quantities. 

All the relevant simulation parameters are summarized in Table 121 Among 
these parameters, the spatial lattice spacing as in physical units correspond- 
ing each /? has been obtained in Ref. [26], together with the corresponding 
parameter u^] the parameter Vc and Vud has been obtained in Ref.[23j. Fi- 
nally, the largest hopping parameter, K^max is chosen such that no exceptional 
gauge field configurations are encountered. This corresponds to the lightest 
pion mass of about 500-600MeV in our simulation. Note that our calculation 
is performed using three set of lattices whose physical volume are about the 
same but with different lattice spacings. The physical size L in the shorter 
spatial direction is about 1.6fm and, for the lightest pion mass in our 
simulation, this gives vtit^L ~ 5 and therefore finite volume corrections which 
might spoil the validity of Liischer's formulae are expected to be small. Three 
different lattice spacings allow us to extrapolate our final results to the con- 
tinuum limit. In order to find the physical point for the charm quark and to 
facilitate chiral extrapolation, six nearby values are taken for both k"'^ and k'^ 
around k^^^ and K^nax^ respectively. 



^.^ Checking single-particle spectrum and dispersion relations 

After inserting a complete set of states, any single-particle correlation function 
can be written in the following form: 

C{t) = (O(t)Ot(O)) = Cne-"""' (24) 

n 
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We may define the effective mass function M^s. (t) as follows: 



M=»,W = log^. (25) 

which in the large temporal limit is dominated by a constant which is the 
mass of the lowest energy gap. Therefore, fitting the effective mass function 
to a constant in a plateau region yields the lowest energy gap. 

In this study, we have calculated single-particle correlation functions of several 
mesons: D*, Di, 7]^, J/ip and tt. D* and Di are the main objects that we want 
to study, other particles are for heavy quark mass interpolation and light quark 
mass extrapolation (chiral extrapolation). As we explained in the previous sec- 
tion, single particle correlation functions for D*, Di, rjc and J/\l/ are measured 
for both zero and non-zero three- moment a. For the pion, only zero-momentum 
correlation is measured. Since we have used anisotropic lattices with enhanced 
temporal resolutions, we obtain descent plateaus for single-particle correlation 
functions. 



Let us first examine the mass of the D* and Di mesons. After obtaining the 
mass values for them under various quark mass parameters ) , the mass 

of the r]c and J/\l/ are used to fix the physical charm quark hopping parameter 
n'^. For this purpose, we demand that the combination ^m^^ + ^mj/^ (spin- 
averaged charmonium mass) reproduces its physical value with the scale set 
by the lattice spacing. This procedure is shown in Fig. [1] where the mass values 
of D* and Di are shown as functions of the charm quark mass parameter k'^. 
The left panel is for D*, and the right one is for Di. In each figure, the 6 
data points with both x and y error-bars are original data for the D meson 
mass. The interpolated point with only y error-bar is the result of charm 
quark mass interpolation. Such interpolations are performed for each light 
quark mass parameter k^'^ although in Fig. [1] one particular light quark mass 
parameter k*"^ is shown. 

After heavy quark mass interpolation for each light quark mass parameter k"*^, 
the mass of the D* and Di mesons {ttid* and moi, respectively) are extrapo- 
lated versus towards the chiral limit = 0. Since our simulation points 
are still far from the true chiral region, we adopted either linear or quadratic 
functions in according to the behavior of the data. This procedure is shown 
in Fig. H 

After all these interpolations and extrapolations, we obtain the mass of D* 
and Di for each lattice spacing. Finally, a continuum extrapolation can be 
carried out for rriD* and m^)^ with linear function in since we are using an 
improved fermion action. This is illustrated in Fig. [31 Our final results for the 
mass of the D mesons are: 
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mo. = 2.008 ± 0.039GeV, mo, = 2.422 ± 0.024GeV, 



(26) 



from which we can see that our results are compatible with experimental result 
within error bars. 

Let us now move on to the dispersion relations for D*, Di, r]c and J/\E'. The 
aim for this study is to verify that we do get single-particle asymptotic states. 
For this we need to know the energy of these particles with definite three- 
momentum. Since the correlators with non-zero three- momentum C (t, k) is 
much noisier than the one with zero three- momentum C{t,0), it is difficult 
to obtain the plateau of the energy directly, particularly for the axial-vector 
meson Di. To get around this, we form the following ratio: 



where 6E{k.) = E(k) — E{0) designates the "kinetic energy" of the particle. 
It turns out that, by forming this ratio, most of the noise is suppressed and a 
plateau for 6E{]<.) can be extracted from: 



These plateaus are illustrated in Fig. HI In fact, for the D*, J/ip and rjc mesons, 
we can also get the plateau directly. There is no need to form the ratio. How- 
ever, if we do form the ratio, the results we get from this ratio are fully com- 
patible with what we get by direct extraction from the original correlators. 
For the axial- vector meson Di, however, forming the ratio helps to suppress 
the noise and to develop the mass plateau. 

After getting the results of 6E{k), the results for -E'(k) = E{0) + 5E{k) is also 
obtained from which one can check the dispersion relation at low-momenta: 

E\k) = m^ + Z-k'' + --- , (29) 



where Z is a parameter to be fitted. In order to recover usual continuum dis- 
persion relation with Z = 1, one has to tune the bare speed of light parameter 
u in the fermion action. This has been done in Ref. for several values of 
f3. In this study, we use the results of u in Ref. [2^ as input parameters. Since 
we have get the value of E(k) in our calculation, we can fit the our data using 
Eq. (l29!l . and get the value of Z. As an illustration, the results of dispersion 
relations are shown in Fig. [5] for certain input quark mass parameters. The 
values of Z are also indicated in each panel. From these results, it is seen that 
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the values of Z is approximately equal to 1 for r^c and J/\&, which suggests 
that our choice for the value of v is approximately right. 



After all these checking, including the mass spectrum and the dispersion re- 
lations for the D* and mesons, we are confident that our finite box can 
accommodate well-established single-particle asymptotic states and we may 
now proceed to study the scattering of D* and D\ mesons at low-momenta. 

^.5 Results for the scattering length and effective range 

As we argued in Sec. [3l only one correlation matrix C{t) have to be computed 
for each symmetry channel of the two-particle system. For each symmetry 
channel, we have studied 5 different non-zero momentum modes. Therefore, 
including the zero-momentum mode, for each symmetry channel, the correla- 
tion matrix for the two-particle system is a 6 x 6 matrix. 

To extract the two-particle energy eigenvalues, we adopt the usual Liischer — Wolff 
method [13]. For this purpose, a new matrix Q(t,tQ) is defined as: 



where to is a reference time-slice. Normally one picks a to such that the signal 
is good and stable. The energy eigenvalues for the two-particle system are 
then obtained by diagonalizing the matrix Q{t,to)- The i-th eigenvalue of the 
matrix has the following behavior in the large (t — to) limit: 



Therefore, the exact energy Ei can be extracted from the effective mass plateau 
of the eigenvalue Aj. 

The real signal for the eigenvalue in our simulation turns out to be so noisy 
that reliable plateau cannot be found directly. Therefore, the following ratio 
was attempted: 



where CD*{t) and Co^it) are one-particle correlation function with zero mo- 
mentum for the corresponding mesons. Therefore, 5Ei is the difference of the 
two-particle energy with the threshold of the two mesons: 




(30) 




(31) 




(32) 
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5Ei = Ei- rriD* - 



(33) 



By taking this ratio, the signal to noise ratio is greatly enhanced. The energy 
difference 5Ei can be extracted reliably from the following effective mass: 

For the A^P channel with correlation matrix given in Eq. (!20|) . the situation is 
illustrated in Fig. Ofor lattices aX (3 = 2.5, 2.8 and 3.2. Six different plateaus 
in each panel correspond to different modes and this procedure is carried out 
for each pair of quark mass parameters (k*"^, k'^). 

With the energy difference 5Ei extracted from the simulation data, one utilizes 
the definition: 



m^* +k2 + y'mf)^ +k2 = + mD* + (35) 

to solve for k^ which is then plugged into the modified Liischer's formula (i.e. 
Eq. ([7]). Close to the scattering threshold, the quantity fc/tan5(fc) has the 
following expansion: 



^ ^ + ^r^e + ■ ■ ■ , (36) 



tan5(A;) 2 

where is the scattering length and tq is the effective range. The l.h.s of 
Eq. (!36l) can also be calculated using Liischer's formula. Therefore, we can fit 
our data with Eq. from which the values of and tq are obtained. Since 
Eq. (15^ is only valid when k is small, we use the data for the lowest 4 modes in 
the fitting. For a particular choice of ), this fitting procedure is shown 

in Fig. [7] for three values of (3 in our simulation. 

After getting the value of and tq for each pair of quark mass parameter 
the results are interpolated versus to the physical charm quark 
mass which is determined by the experimental value of ^m.^^ + |mj/^. This 
is shown in Fig. [HI The interpolated data are then taken for the chiral ex- 
trapolation. In this step, the results for and rg are extrapolated versus 
towards the chiral limit as shown in Fig. [91 Finally, continuum limit is taken 
by a linear extrapolation in for the results of Cq and Tq obtained after chiral 
extrapolation. The final results for the scattering length oq and the effective 
range tq in this channel is shown in figure [lOl After these extrapolations, we 
obtain the scattering length qq and the effective range Tq for the Ai channel: 



Co = 2.53 ± 0.47fm, tq = 0.70 ± O.lOfm . (37) 
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2.5 


-0.026(0.003) 


5.23(0.65) 


2.8 


-0.064(0.005) 


0.16(0.18) 


3.2 


-0.053(0.016) 


0.92(0.93) 



Table 3 

Results for the lowest and the corresponding values for cot a{q) as given by Eq. ([8]) 
for different values of /3 in the simulation. Corresponding errors for the quantities 
are also given in the parenthesis. 

This result is for A^^^ channel which, in the notion of continuum quantum 
numbers, corresponds to = 0^. The signal in other channels is much noisier 
than that of Ai channel and it seems that more statistics and/or better 
interpolation operators are needed for a reliable extraction of the scattering 
parameters. Results for the scattering length oq and the effective range tq at 
various light quark mass parameters for three values of (3 are also listed in 
Table m for reference. The results after the chiral extrapolations and the final 
results in the continuum limit are also shown in the table. 



4-4 Possibility of a shallow bound state 

To explore the possibility of a bound state, we recall that for a bound state 
to exist, has to be negative and in fact q"^ —>■ —oo a.s L oo. This results 
in the condition: cot(j(g) = —1 as discussed in the subsection 12.11 Eq. On 
the other hand, a scattering state will have: ~ (1/L) as L ^ oo. Results for 
the lowest (negative) and the corresponding values of cot a{q) as computed 
from Eq. ([8]) are listed in Table [31 It is seen that our results for cot cr(g) for the 
lowest (negative) are all positive. The absolute values for the lowest are 
also not large. Our results obtained so far seems to be more consistent with a 
scattering state than a bound state. 

One could investigate this possibility from another point of view, namely by 
the values of scattering length and effective range. The value of effective range 
ro obtained is much less than the size of our box so that using Liischer's for- 
malism is justified. Since we are studying the scattering near the threshold, it 
is appropriate to study the problem using non-relativistic quantum mechanics. 
Within non-relativistic quantum mechanics, it is known that, if a shallow 
bound state emerges in s-wave potential scattering at low-energies, the scat- 
tering length of the system will diverge. In fact, if the potential acquires an 



^ See, for example, "Quantum Mechanics (non-relativistic theory)", 3rd ed., 
L.D. Landau and E.M. Lifshitz, Pergamon Press, §133. Note also that our defi- 
nition on the scattering length differs from theirs by a sign. 
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infinitely shallow bound state, the scattering length should approach negative 
infinity |l6j. Our lattice results for the scattering lengths indicate that it is 
quite large but positive. This usually happens when the potential is on the 
verge of developing a shallow bound state. Note that this argument is generally 
valid for a wide variety of potentials. 

If we further approximate the potential by a square-well potential, we could 
even estimate the depth Vq and the range of the potential R from our lattice 
results on oq and tq. We find that, R = tq = 0.70(10)fm and Vq = 73(21)MeV. 
These values for a square-well potential also gives no bound states. If we fix 
tqR = 0.7fm, the first bound state will occur at about Vq ~ 92MeV. 



5 Conclusions 



In this paper, we present our quenched anisotropic lattice study for the scat- 
tering of D* and Di mesons near the threshold. The calculation is based on a 
finite-size technique due to Liischer which enables us to extract the scattering 
phases from the exact two-particle energies measured in Monte Carlo simu- 
lations. Our study focuses on the s-wave scattering in the channel = 0~ 
and the scattering threshold parameters, i.e. scattering length and effec- 
tive range tq are obtained. After the chiral and continuum extrapolations, we 
obtain: oq = 2.53(47)fm and tq = 0.70(10)fm, indicating that the interac- 
tion between a D* and a Di meson is attractive in this channel. As for the 
other channels, although we have also computed the correlation matrices, but 
the signal is too noisy to obtain definite results. Better operators and more 
statistics are probably needed in further studies. 

Based on our results for the scattering phases near the threshold, we have also 
discussed the possibility of a shallow bound state in this channel. We inves- 
tigate the quantity cot cr which should approach (—1) for a bound state. Our 
results for this quantity are all positive. Our results for scattering length are 
also positive. Based on these indications, it seems that, although the interac- 
tion between the two charmed mesons is attractive, it is unlikely that they 
form a genuine bound state right below the threshold. The lowest two-particle 
state is likely to be a scattering state. This result might shed some light on the 
nature of the recently discovered Z+(4430) state by Belle. However, we should 
emphasize that, our lattice calculation is done in a particular channel only 
and it is within the quenched approximation. Obviously, to further clarify the 
nature of the structure Z"'"(4430), lattice studies in other symmetry channels 
and preferably with dynamical fermions are much welcomed. 
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Table 4 

Results for the scattering length ao and the effective range vq at various light quark 
mass parameters for three values of fi. The results after the chiral extrapolations 
and the final results in the continuum limit are also shown. 
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